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The intuitiveness of the tensor network graphical language is becoming well known 
through its use in numerical simulations using methods from tensor network algo- 
rithms. Recent times have also seen rapid progress in developing equations of motion 
to predict the time evolution of quantum entanglement [Nature Physics, 4(4) :99, 
2008] . Here we cast these recent results into a tensor network framework and in do- 
ing so, construct a theory which exposes the topological equivalence of the evolution 
of a family of entanglement monotones in arbitrary dimensions. This unification was 
accomplished by tailoring a form of channel state duality through the interpretation 
of graphical tensor network rewrite rules. The introduction of tensor network meth- 
ods to the theory of entanglement evolution opens the door to apply methods from 
the rapidly evolving area known as tensor network states. 
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There has been increasing recent interest in the topic of evolution of quantum 
entanglement — inspired by the work Konrad et al. considered the evolution 
equation for quantum entanglement for qubits [l[ and its later extension to higher 
dimensions can be found in @, Q. The multipartite case was characterized by 
Gour \4\. Before people have generally dealt with entanglement evolution 
on a case- by-case basis. In this spirit, several authors considered non-Markovian 
entanglement evolution Similarly, upper and lower bounds can be obtained 

for special classes of states such as graph states 

Our approach aims to develop a tensor network diagrammatic method and to 
show which aspects of entanglement evolution can be described in those terms. We 
found that key results of Tiersch [2[ can be reproduced in the language of tensor 
networks with several appealing features. Konrad et al. themselves developed and 
relied on a sequence of diagrams, that aided in their exposition of the topic at 
hand [l|. These diagrams were utilitarian in nature and used to convey a point. 
However, there exists a diagrammatic graphical language, of which Hilbert spaces 



are a complete model, known as Penrose tensor networks or string diagrams [10 



These diagrams can aid in intuition and also represent mathematical equations 
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The diagrammatic language built on to describe the tensor networks dates back 



to Penrose [10J. Recently, we have also done work on the theory and expressive- 
ness of tensor network states I17H . What we add to the literature on tensor 
network states is the development of some theory which connects these ideas to a 
specific version of channel-state duality. This offered an advancement and exten- 
sion of the prior graphical aid of [l| • Through this rigorous graphical language a 
unified common topological structure of the networks representing the evolution 
of entanglement also emerges. 

The introduction of a tensor network theory of entanglement evolution opens 
the door to apply numerical methods from the rapidly evolving area known as 



tensor network states [151 . 1161 . Il8| . In recent times, tensor networks states have 
been successfully used to simulate a wide class of quantum systems using classical 
computer algorithms. In addition, theory has been developed to use these methods 
to simulate classical stochastic systems [l8j]. At the heart of these methods are 
tensor network contraction algorithms, pairing physical insight with a new theory 
behind making approximations to truncate and hence minimize the classical data 
needed to represent a quantum state. The key properties and functions of these 
algorithms are known to be described in the graphical language we have adapted 



to entanglement evolution [lOj, [12|, [13 



Structure. We will begin in Section |TT] by reviewing those tensor network com- 
ponents that we will use. This introduction paves the way for a review of the so 
called snake equation. This equation is then used to prove a graphical version of 
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map-state duality — inspired by and tailored specifically to the problem at hand. 
Section II III then applies these building blocks and associated tools from tensor 
network states to the theory of entanglement evolution. We begin by recalling 
the fundamental definitions and a quick outline of the known theory, we then cast 
these methods into our framework and pinpoint a wide class of measures which all 
share the same network topology in the tensor network language. The conclusion 
presents an overview of potential future directions by listing some open problems 
and new directions. 



II. PENROSE GRAPHICAL NOTATION AND MAP-STATE DUALITY 



We begin by showing how quantum states and operators may be represented 
graphically using the tensor network language. We then proceed to more complex 
topics, including the state-map duality and its graphical equivalents. 

It is typical in quantum physics to fix the computational basis and consider a 
tensor as an indexed multi-array of numbers. For instance, 

(a) (b) I* 



represent the tensor (a) ipi and (b) T l - k . The scenario readily extends. 

Graphically a two party quantum state would have one wire for each index. 
We can write in abstract index notation as tpij or in terms of the Dirac notation 
convention (which we mainly adopt here) as \ifj) = Yliij c ij K) li)- Appropriately 
joining a flipped (transpose) and conjugated (star or overbar) copy of the quantum 
state allows one to represent the density operator p = c ijC kl (kl\ as 
follows 0. 




FIG. 1: Penrose's graphical representation of a pure density matrix. 



Here a box with two inputs (i, j) and two outputs (A;, /) would be used to rep- 
resent a general density operator p = ^ P%j kl with no known additional 
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structure. Those familiar with quantum circuits will recognize this graphical lan- 



rcm 



guage as their generalization [13| 



Remark II. 1 (Tensor valence) A tensor with n indices up and m down is called 
a valence-(n, m) tensor and sometimes a valence-k tensor for k = n + m. 

Remark II. 2 (Diagram convention — top to bottom, or right to left) 

Open wires pointing towards the top of the page, correspond to upper indices 
(bras), open wires pointing towards the bottom of a page correspond to lower in- 
dices (kets). For ease of presentation, we will often rotate this convention 90 
degrees counterclockwise. 

We are interested in the following operations: (i) tensor index contraction by 
connecting legs of different tensors, (ii) raising and lowering indices by bending a 
leg upwards or downwards, respectively and thus taking the appropriate conjugate- 
transpose and (iii) a duality between maps, states and linear maps in general, called 
Penrose wire bending duality. This duality will play a major role in our application 
of the language to entanglement evolution. The duality is obtained by noticing 
that states are tensors with all open wires pointing in the same direction. Bending 
some of the wires in the opposite direction makes various linear maps out of a given 
state. Given the computational basis, this amounts to turning all kets belonging 
to one of the Hilbert spaces into bras and vice versa. A bipartite state in the fixed 
standard basis 



1,3 



for example, is dual to the operator 

M = £A>|i>01- (2) 



h3 



This result can be shown by considering components of the equation (A ® 1) |$ + ) 
where |$ + ) is the Bell state ($5§. 

The Choi-Jamiolkowski isomorphism can be cast of Penrose Duality. 

To do this we must recall Penrose's cups, caps and identity wires. As in [To| , these 
three tensors are given diagrammatically as 

(a) & (b) □ * (c) n * 
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By thinking of these tensors now in terms of components, e.g. Sy is 1 whenever 
i = j and otherwise, we note that 

i = w 01 = E !*><*! (3) 

ij k 

- (11| + ■ ■ ■ + (nn\ = (til = E ( kk \ ( 4 ) 



ij k 

|00) + |11) + ■ ■ ■ + \nn) = J2 5 H \ti) = E \ kk ) ^ 

ij k 

where the identity map (a) corresponds to Equation (j3J), the cup (b) to (TjJ and 
the cap (c) to (j3J). The relation between these three equations is again given by 
Penrose's wire bending duality: in a basis, bending a wire corresponds to changing 
a bra to a ket, and vise versa, allowing one to translate between Equations ([3]), 
flU) and (|5|) at will. 

The contraction of two tensor indices diagrammatically amounts to appropri- 
ately joining open wires. Given tensors T l j k , A 1 n and B m we form a contraction 

by multiplying by Sj 5 k resulting in the tensor 

T l jk^ 3 iB k m :=T\ m (6) 

where we use the Einstein summation convention (repeated indices are summed 
over) and the tensor T\ m is introduced per definition to simplify notation. In 
quantum physics notation, this is typically expressed in equational form as 

r = rSm \lm) (*| = E T l jk A\B k m \lm) (*| . (7) 

ilm ijklm 



Result II. 3 (Graphical trace — Penrose [10]) Graphically the trace is per- 
formed by appropriately joining wires. The following depiction in Figured (b) is 
Penrose 's representation of a reduced density operator, where the stars on the ^> 's 
represents complex conjugation. 



We have now presented the key tensor network building blocks used in this 
study. In practice, tensor networks contain an increasing number of tensors, mak- 
ing it difficult to form expressions using (inherently one- dimensional) equations. 
The two-dimensional diagrammatic depiction of tensor networks can simplify such 
expressions and often reduce calculations. Here we tailor this graphical tool to 
depict internal structure and lend insight into presenting a unified theory of en- 
tanglement evolution. A key component of this unified view relies on the natural 
equivalence induced by the so called snake equations, which we will review next. 
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Tr 2 (p) = ^ 



I w 



FIG. 2: Graphical representation of a partial trace. The network (b) is due to Penrose. 
(Note that Penrose used reflection across the page to represent adjoint, yet we have 
placed a star on \& to represent conjugate for illustrative purposes.) 



The snake equation [10J amounts to considering the transformation of raising 
or lowering an index. One can raise an index and then lower this index or vice 
versa, which amounts essentially to the net effect of doing nothing at all. This is 
captured diagrammatically by the so called, snake or zig-zag equation, as 



c 



together with its vertical reflection across the page. In tensor index notation, it is 
given by 5ji5 lh = 8-*. 

As a tool to aid in our analysis of entanglement evolution, we have noticed a 
succinct form of map state duality that works by considering the natural equiv- 
alence found by the snake equations inside an already connected diagram. We 
expect that this formulation and its interpretation will have applications outside 
of the area of entanglement evolution. To date, the snake equation has already 



been used to model quantum teleportation [19|, [20| and so the interpretation here 
is a different one. 

Lemma II. 4 (Graphical map-state duality) Every bipartite quantum state 
is naturally equivalent to a single sided map [ip] acting on the Bell state |$ + ) 
with graphical depiction as follows. 

We will make a few remarks related to Lemma 111.41 To prove the result using 
equations (instead of diagrams which was already done) we proceed as follows. To 
consider \ifj) as a map [if}], we act on it with a Bell state as [ip] := (1 (g) ($ + |) \ip). 
(Notice that the net effect is to turn kets of one of the spaces into bras and vice 
versa, i.e. |0) <g> |0) becomes |0) (0|.) We then note that \ip) = ([ip] <g> 1) |$ + ) and 
hence that \ip) = ([(1 Cg> ($ + |) \ip)] ® 1) |$ + )- The result seems more intuitive by 
considering the graphical rewrites. 

To explain the diagrammatic manipulations, begin with Figure E] (a), represent- 
ing the state \ip). We can think of this (vacuously, it would seem) as a state being 
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FIG. 3: Graphical channel state duality of a state acted on by the identity channel, (a) 
The state can be thought of as being acted on by the trivial or identity channel on 
both sides, (b) From the snake equation an identity wire can be equivalently replaced 
with its 'S' shaped deformation. This can then be thought of as a map [ip] (enclosed 
by a box with dashed lines on its edges) acting on a Bell state (enclosed with a triangle 
with dashed lines on its edges), (c) An equivalent expression found by redrawing the 
shape of the tensor representing \ip), into a box. 



acted on by the identity operation. Application of the snake-equation to one of 
the outgoing wires allows one to transform this into the equivalent diagram (b). 
We can think of (b) as a Bell state (left) being acted on by a map found from 
coefficients of the state We have illustrated this map acting on the bell state 
in (b) by a light dashed line around \ip). We then rewrite the figure, capturing the 
duality and arriving at (c). 



Remark II. 5 (Properties of the map if)) The map [ip] := (1 Cg> ($ + |) evi- 
dently has close properties with the quantum state \ip). For instance, Trff//] = 
(ip\ifj). The state \ip) is invariant under the permutation group on two elements 
iff W\ = H>] T ■ The state is SLOCC equivalent to the generalized Bell state 
iff det[ifj] and the operator rank of [ifi] gives the Schmidt rank of Rank 
is a natural number and although it is invariant under action of the local unitary 
group, it is not what is called a polynomial invariant or algebraic invariant since 
rank is not expressed in terms of a polynomial in the coefficients of the state. The 
eigenvalues of the map [ip] are all that is needed to express any function of the 
entanglement of the bipartite state 



Remark II. 6 (Extending Lemma |II.4j) One can readily extend Lemma II. 4 
and consider instead the scenario given in the following Figure. 



Ill ENTANGLEMENT EVOLUTION 






FIG. 4: Graphical formulation of channel state duality, (a) The diagram representing 
($ ® 1) (b) The tensor network of the same state and the application of the natural 
equivalence induced by the snake equation per Lemma III. 41 (c) We then think of the 
state tp as a map Yl^ij K) 01 ac ting on the state |$). 



III. ENTANGLEMENT EVOLUTION 



Remark III.l (Background reading) For more details on standard results in 



the area of open systems theory, see for instan ce 12 Ah For result results unifying 



several pictures of open systems evolutions see Xl Ti l. 



Quantum operations are maps acting on quantum states p such that they pre- 
serve the positivity and boundedness of the corresponding density operator [2l|. 
Such maps are known as completely positive. They can be used to represent the 
action of unitary maps or, more generally, decoherence and measurement or com- 
binations thereof. In general a completely positive map $, called a superoperator, 
may be described in terms of Kraus operators A k , 



$[p] =Y,A kP A 



The map $ is said to be trace-preserving iff Tr(p) = Tr($p) for all positive operators 
p and it can be shown that this is true iff Ylk ^-fc^-l = 1- Such maps faithfully 
preserve the normalization of p. 

We consider bipartite states p G Ha^Hb- Since such states may be entangled, 
one is often particularly concerned with how entanglement evolves or changes under 
completely positive maps. Thus, given an entanglement measure C and a local 
completely positive map $ = $^ <g> § B acting as $ [p] = ^ k (A k <g> B k )p(A\ <g> B\) 
we will want to compute C[($,4 ® $b)p], for some entanglement measure C. An 
elementary, but still critically important, example of such maps are one-sided maps 
of the form $ = $^ ® 1b- 



S 
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Entanglement evolution under one-sided maps has been characterised in the 
literature. In [l| Konrad et al. showed that for pure bipartite qubit states 

c[{*a ® i) m w] = c [$a ® i)i<f + )(<f + i] -c[m w] (9) 

and so they factored the function C[$, l^)] of two variable arguments into a product 
of functions, one depending only on the starting state and the other only on the 
completely positive map characterising the one-sided evolution. Thus, calculating 
the evolution of entanglement for a single initial state, |$ + ), allows us to compute 
it for all other pure initial states. As will be seen later, this provides an upper 
bound on the quantity of evolved entanglement for all initial mixed states. Here 

|$+) = -^(|00> + |ll» (10) 

is again the typical Bell state and C is the concurrence. For mixed states they 
found the upper bound to be given by the product [lj 

C[($ A <8>l)p] <C[($ A ®l)|$ + >($+|] -C[p\. (11) 

In [2| the analogous result was obtained for higher- dimensional systems, but 
where C was replaced by the (/-concurrence, defined for (i-dimensional pure states 

d 

|V> = Yl A * N) ® \3) (12) 

with A := [ip] as 

C d U)} = d- [det(AU)] 1/d (13) 



For mixed states the ^-concurrence is defined through the convex roof construc- 



tion 22 



g d [p] = inf ^2piC d [\ipi)] 



(14) 



where the infimum is taken over all convex decompositions of p into pure states, 
i.e., all decompositions of the form 



P 



\A) (fa 



(15) 



with p.i > 0, YliPi = 1 an< i the need not be orthogonal. Hence, we note in 
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particular that it is not sufficient to consider only the spectral decomposition. 

The results on factorisation of entanglement evolution were further extended to 
the multipartite case in Q , where Gour showed that the same factorization result 
holds for all SL-invariant entanglement measures. These are defined as: 

Definition III. 2 (SL-invariant measures |3. |2~2|) Let G be the group G = 
SL(di, C) (g> SL(d 2 , C) <E> ■ ■ ■ <S>SL(d n , C). The groups SL((4, C) are the special linear 
groups represented by alldkxdk matrices with determinant 1. A SL-invariant mea- 
sure of multipartite entanglement C : B(H) — > M + maps the bounded operators on 
the Hilbert space % (not necessarily normalised density matrices) to non-negative 
real numbers. C must satisfy the following properties: 

1. For every g G G and p G B(Ji) we have that C[gpg^] = C[p\. 

2. It is homogenous of degree 1, C[rp] = rC[p]. 

3. For all valid mixed states it is given in terms of the convex roof construction: 

C[p]=mmJ2PiC[\ipi)m]. (16) 

i 

Here the minimum is taken over all convex decompositions of p, i.e. all 
decompositions into p = YliPi 1^) i^il where (ipi \ ipj) is not necessarily zero. 

Now we will show how the above results, both those for mixed and pure states 
and generally for qudits, can be described through the use of Penrose graphical 
calculus and the Choi-Jamiolkowski isomorphism. More precisely, we will show 
that 

C[($ A ® 1) m (VI] = C[($ A ® 1) |$ + > <$ + |] ■ CU) M]. (17) 

The factorised form of the entanglement evolution is desirable as it simplifies and 
unifies calculations of entanglement once a one-sided map has been applied. To 
arrive at the factorisation formula we first require one further result. 

Given annxn matrix A with det(A) ^ and an SL-invariant n-dimensional 
entanglement measure C, we find that 

C[ApA^] = C[\ det{A)\ 2 ' n gpg ] ] = |det(A) 2 / n | C[p\. (18) 

where we have defined A := det(A) 1 ^ n g, with g G SL(n). 

As will be seen, this calculation (I18p will be applied to maps A that arise under 
wire bending duality from states. Consider 
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and by Lemma III. 41 we establish equality of (a) and (b). Hence, every bipartite 
quantum state 1^) can be represented as the generalized Bell state |$ + ) acted on by 
a single sided map [ip]. Let A := [ifj] then the rank of A determines entanglement 
properties of A is full rank iff is SLOCC equivalent to |$ + ). In fact, A is 
full rank iff det(A) 7^ or equivalently iff det(A^A) > 0. This will soon become 
relevant as we will use this duality to view \if>) as a single sided map and rely on 
the determinant to evaluate the entanglement measure of interest. 

Now that we have established a factorisation of the above formula (ITS]) we will 
show that the factor det(A) 2//n is nothing other than C[($ <8> 1b) \4> + ) Let us 

then consider evolution of the bipartite state under the superoperator $, with 
a tensor network given as 




We then apply Lemma 111.41 to both sides and arrive at 




So \ip) becomes a superoperator acting on the density operator 



k 



\A k ) (A fe | 



(19) 



which can be equivalently expressed in Figure |5](c) where the remaining subfigures 
(a, b) summarize the duality behind entanglement evolution. Hence, Lemma III. 41 
allows us to show the equivalence of the states (a) and (c) in Figure H] and from 
this it follows that 

($<g)l)|V><V| = (l®$v)P$- (20) 

which is alternatively expressed as (1 (g) [if)]) J2k l-^-fc) (-^fcl- 

Furthermore, notice that any state may be written using the Schmidt de- 
composition as = J2 k ^JuJk \kk). Thus, the map may be considered simply 
as an action of a matrix 



V^k\ k ) ( k \ 



(21) 



where the factor of y/n enters because of normalization. Because of the purity of 
j^) we arrive at only a single Kraus operator. Thus we find that 



C[($®1) |V) {M\ = det(A) 2 /"C[p $ ] 



(22) 
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Notice, that det(A) 2 /™ = nY\^ =1 u\/ n . This is exactly the definition of Q- 
concurrence, reflecting the fact that ^-concurrence is the only SL-invariant measure 
of entanglement for bipartite qudits. The second factor in the above equation (122"]) 
is simply equivalent to 

C[p $ ]=C[($®l)|$ + }($ + |] (23) 

where 

1 n 

l $+ > = 7^El^>- (24) 

V k=l 

We can therefore rewrite the equation as 

C[($ ® 1) M = C[($ ® l)|$ + )($ + |] • C[|^> (VI] • (25) 

n 

= C[($®l)|$ + >($+|].nn^ (26) 

fc=i 

The term C[($ £g> l)|$ + )($ + |] still needs to be evaluated. However, its value 
tells us how concurrence evolves under $ ® 1 for all pure states. It thus only needs 
to be computed once. Here the entanglement measure C is now understood to be 
the ^-concurrence. The equation factorises concurrence into a part entirely due to 
evolution dynamics and a part entirely due to the initial state. 

Having established the evolution equation (I25|) . we will now consider an upper 
bound. Since the measure C is given by a convex roof construction for mixed 
states, it is a convex measure — i.e. 

C^PkPk] <5>*C[p*]- (27) 
k k 

We can use the convexity of the measure C to find an upper bound on the entan- 
glement evolution for mixed states. This upper bound is given by 

C[($ ® l)p] < C[($ ® l)|$ + }($ + |] ■ C[p\. (28) 



IV. CONCLUSION 

We have shown how entanglement evolution can be understood and even unified 
in terms of a tailored variant of the Penrose graphical calculus. We developed a key 
graphical tool to simplify the derivation of several important recent result in quan- 
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(a) 



- A - 




- A - 



(b) 



(c) 



- ¥ - 




- ¥ - 



- V - - V - 



FIG. 5: Tensor network summary of the crucial duality behind the theory of entangle- 
ment evolution, (a) State \ip) acted on by a single sided superoperator $. (b) Applica- 
tion of Lemma III. 41 establishes a duality mapping the scenario in (a) into view where 
\ip) becomes a single sided superoperator [ip] which acts on the bipartite mixed state 
p$ = |Afc) (Afc|. (c) An equivalent expression for the valence four tensor p$. 



turn information theory — that of entanglement evolution of bipartite quantum 
states. The result was originally obtained through the Choi-Jamiolkowski isomor- 
phism, which we have shown in the diagrammatic language reduces to bending 
a contracting wire into a snake. We introduced the necessary tensor networks 
and showed how the evolution of ^-concurrence, the only SL-invariant measure of 
bipartite entanglement, can be obtained from these network diagrams. We thus 
developed a tensor network approach that, through the use of Penrose duality, 
recovers the known theory of entanglement evolution. This goal necessitated the 
introduction of several new methods which could prove useful in other areas. In 
particular, we anticipate the existence of other problems that could be simplified 
through the use of the diagrammatic variant of channel state duality we have de- 
veloped for our purposes here. This cross pollination could lead to further insights 
into the strengths and limitations of the approach we have presented here. 

As for the entanglement evolution itself, it is now fairly well understood for 
pure initial states and one-sided channels, few results are known for more general 
maps and for mixed initial states. Progress in this direction through the use of 
the Choi-Jamiolkowski technique is inhibited by the fact that the dual equivalent 
of a mixed state is a super-operator rather than a simple matrix. Mathematically, 
given a mixed state in its spectral decomposition p = YlkPk l"0fe) (V'fcli ^s wire bent 
dual is the map (see part (a) of figure E]) 

% p = Y,PkA k - A[, (29) 

k 

where Ak is the wire bent dual of the state \ipk)- With pure states we have the 
fortune that $ p reduces to a single matrix and we can thus greatly simplify the 
calculation. The fact that it is now a genuine super-operator limits the applicability 
of this approach. One would thus expect that a qualitatively different approach is 
needed to find a closed form expression for entanglement evolution of mixed states. 
And in fact, the very existence of such an expression is not guaranteed. However, 
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finding it would immediately also yield a closed form expression for the evolution 
under maps acting on more than one subsystem $1 <g> $2. To see this, notice that 
$1 <S> $2 = ($1 <S> 1) • (1 ® $2) thus treating the map as a sequence of two maps 
acting on only one subsystem. 

As with the purely algebraic approach, the drawback of the tensor network 
approach is that it relies on the simplifying properties of bipartite SL-invariant 
measures of entanglement. The key contribution being a conceptual aid and expo- 
sure of the underlying mathematical structure at play. Namely, the property that 
acting on a single subsystem of a quantum state with a non-singular matrix as 
(A 1) p (A^ ® l) multiplies an SL-invariant entanglement measure with a func- 
tion of det(A). Other entanglement measures may not satisfy this property, and 
indeed the usual concurrence only satisfies it in two dimensions. 
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